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We investigate basic physical quantities for quenched simulation with domain-wall fermions and the DBW2 
gauge action. Masses and decay constant of pseudoscalar mesons are measured. Scaling properties are tested. 



1. INTRODUCTION 

Domain wall fcrmion (DWF) |^J^ is a promising 
tool to calculate hadronic quantities where chi- 
ral symmetry is important. Precise calculation of 
kaon mixing parameter Bk is one success of this 
formulation 1^. More challenging calculations |^,|| 
of kaon weak matrix element for e'/e show a first 
reliable determination within lowest order chiral 
perturbation theory. Attempt to calculate nu- 
cleon matrix elements is also underway]^. 

The good chiral property of DWF comes at an 
obvious cost, which is a substantial increase of 
computational effort associated with the size of 
the fifth dimension L^. This makes taking the 
continuum limit by performing lattice finer than 
a"^ ~ 2 GeV difficult. Going to a coarser lattice 
to have another point for the continuum extrap- 
olation could be a better strategy. However, re- 
sults on coarser lattices show that an incredi- 
bly large size for the fifth dimension is required to 
keep the residual chiral symmetry breaking small. 

Motivated by the fact that improving the 
gauge action reduces the residual chiral symme- 
try breaking at finite Lg significantly at a"^ ~ 2 
GeV0, we found the DBW2 action|||l, which is a 
non-perturbatively determined Iwasaki-type ac- 
tion[^o|, gave much more improvement for the 
chiral properties ||ri|,|l^. Here we investigate ba- 
sic physical quantities for the DBW2 gauge action 
for quenched DWF to see the viability of this ap- 
proach down to fl-^ ~ 1.3 GeV. 



Table 1 

Simulation parameters and calibration. 



*Talk presented at Lattice 2001, Berlin, Germany. 
tWe thank RIKEN, Brookhaven National Laboratory and 
the U.S. Department of Energy for providing the facilities 
essential for the completion of this work. 
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1.28(3) 


1.04 
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16 


200 


1.7(1)10-5 


1.97(4) 



We perform the simulation on a 16^ X 32 lat- 
tice with DBW2 action at /3 = 0.87 and 1.04 
in Iwasaki type parameterization, which corre- 
sponds to ~ 1.3 and 2 GeV respectively. The 
domain wall hight M5 is tuned in the first deci- 
mal place so that it minimizes the residual quark 
mass rriresj although it is not necessary that it be 
done to this accuracy. 



2. RESULTS FOR a 



The residual quark mass m 
through the ratio 

R{t) = {J5,{t)p{t)) I {p{t)pm 



1.3 GeV 

res is defined 



(1) 



at t ^ a, where Jag represents explicit breaking 
of chiral symmetry in the the divergence of "con- 
served current" A^, and P is the pseudoscalar 
density operator. We list the value of rrires in Ta- 
ble |l|. For fl-^ ~ 1.3 GeV unrenormalized value 
of nires — 0.7 MeV, which is smaller than that 
obtained for the Wilson gauge action at finer lat- 
tice, a~-^ ~ 2 GeV and the same Lg, rrires — 2.4 
MeV§. 

The pion mass can be extracted from a two 
point function constructed from pion interpolat- 
ing operators. We consider two different two- 
point functions. One is the pseudoscalar - pseu- 
doscalar, (PP), the other is the local axial - ax- 
ial correlator, {AoAq). The former has a 1/m^ 
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Figure 1. Pion mass squared as a function of 
quark mass for a^^ ~ 1.3 GeV. Asterisk indicates 
the point ruf = —rrires- 

and l/m pole from zero modes in finite volume, 
while the later has only l/m. Figure ^ shows 
the pion masses from those different two point 
functions are consistent with each other down to 
nif = 0.02. The difference begins to show up at 
the lightest mass mfa — 0.01, however it is quite 
small compared to the value of mass itself. 

A linear fit of the (AqAq) correlator for 0.01 < 
m/ < 0.06 is quite good {x^/dof = 0.03/3). 
However, x-intercept is m/ = —1.3(4)10"'^, so it 
overshoots to/ = —rrires, where the pion mass 
should vanish. This location of the intercept is 
guaranteed with accuracy of 0{a^) by the Ward- 
Takahashi identity when all possible finite size ef- 
fects have been removedl^. With this constraint 
one can fit to a quenched chiral log formulae, 

(totto)^ = aoTO (1 + ai logTO) + a2TO^ (2) 
m = m/a + TOresO. (3) 

Results are listed in the Table |[ Both fits with 
or without quadratic term give reasonable value 
of ai. 

Table 2 

Fit result of pion mass. "-" indicates that the 
parameter is constrained to be zero. 
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Figure 2. Pseudoscalar decay constant as a func- 
tion of quark mass at ~ 1.3 GeV. 

We determine the pseudoscalar decay constant 
using two methods [| I . One is an estimate from 
the {AqAq) correlator with the renormalization 
factor Za nonperturbatively determined from the 
ratio of (^o^*) and {AqP). The other uses 
the {PP) correlator with the equation A^^^ ~ 
2 (to/ -I- mres)P from the Ward-Takahashi iden- 
tity. As seen in the Figure |2| the two methods give 
consistent results, indicating the effect of chiral 
symmetry breaking can be absorbed in the shift 
of TO/. Also the possible zero mode effects, which 
influence the two correlators differently, appear 
to be small, which is a confirmation of what we 
saw for the pion mass down to to/ = 0.02. 

3. SCALING PROPERTIES 

Figure || plots the rho meson mass extrapolated 
to the chiral limit as a function of lattice spac- 
ing. Both axes are normalized by rp calculated 
from the heavy quark potential. Together with 
the data of this work, results with Wilson gauge 
actional are shown. Those data show that the 
scales determined by different methods agree very 
well for both the DBW2 and the Wilson gauge ac- 
tion for DWF. This good scaling is also seen in 
the nucleon mass and K* mass normalized by rho 
meson mass. This also holds when one chooses rg 
for the input of scale, which is obvious from above 
discussion. 

Scaling of the pseudoscalar decay constant 
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Figure 3. Consistency check of scales from the 
rho meson mass and the heavy quark potential. 

(Fig. 1^) is also good. is in good agreement 
with the experimental value, fx appears to be 
smaller than the experimental value, which could 
be caused by the effects of quenching |l3|. 

4. CONCLUSION 

We have attempted to use the DBW2 gauge 
action in a quenched domain wall fermion simu- 
lation. Good chiral properties for ~ 1.3 GeV 
has been observed. We have investigated scaling 
of hadronic observables, which also turned out to 
be good. We have also confirmed that the val- 
ues were consistent with those obtained with the 
Wilson gauge action. 

The statistics presented here is not sufficient for 
the precision measurement for those observables 
calculated here. However, we can argue that there 
is no evidence for any difficulty with this choice of 
gauge action, rather it has lots of good features. 
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Figure 4. Scaling plot of nucleon and K* mass. 
Diamonds are experimental values. 
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Figure 5. ScaHng plot of pseudoscalar decay con- 
stant. Diamonds are experimental values. 

An important benefit from the use of this gauge 
action is that residual chiral symmetry break- 
ing is quite small even for a coarse lattice with 
a^^ ~ 1.3 GeV. Therefore it may be suitable for 
demanding calculation like e'/e. Once we have 
the results for ~ 1.3 GeV and finer lattices, 
we will be able to discuss the scaling property for 
those quantities. 
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